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摘 要 : 为 了 研究 液 氧 煤油 发 动机 尾 焰 冲 击 导 流 槽 的 流 场 规律 ， 采 用 CFD 软件 计算 了 液 氧 煤油 发 动机 尾 焰 对 无 导 流 装置 、 
有 镍 形 导 流 装置 和 有 圆锥 体 导 流 装置 三 种 不 同 导 流 槽 的 冲击 流 场 进行 了 三 维 数值 模拟 , 对 比分 析 了 冲击 不 同 导 流 槽 的 尾 焰 燃气 
流动 。 结 果 表 明 : 液 氧 煤油 发 动机 尾 焰 燃气 垂直 冲击 无 导 流 装置 的 导 流 模 ， 沿 着 导 流 槽 侧面 流动 和 向 上 反射 的 燃气 会 对 导 流 模 
底面 正 上 方 环 境 产生 高 温 影 响 ; 狗 形 导 流 装置 能 够 避免 治 着 导 流 槽 侧面 流动 和 向 上 反射 的 燃气 对 导 流 槽 底面 正 上 方 环境 的 高 温 
影响 ;圆锥 体 导 流 装置 能 够 避免 向 上 反射 的 燃气 对 导 流 槽 底面 正 上 方 环 境 的 高 温 影 响 ， 在 避免 燃气 沿 着 导 流 槽 侧面 流动 方面 仍 
有 改进 余地 ;， 相 对 尾 焰 燃 气 冲击 无 导 流 装置 的 导 流 模 ， 狠 形 导 流 装置 和 辆 锥 体 导 流 装 置 表面 最 大 压力 分 别 增 大 了 10. 31% 和 
33. 81%， 负 形 导 流 装置 和 圆锥 体 导 流 装 置 受 尾 焰 燃气 冲击 的 最 高 温度 分 别 降低 了 4. 04% 和 8. 95%。 
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Abstract: Aimed at studying the flow characteristics of LOX/kerosene engine exhaust plume impinging on the 
diversion trough, three-dimension flow fields of LOX/kerosene engine exhaust plume impinging on three diversion 
troughs of without deflector, with a wedge deflector , with a conical deflector were calculated with the CFD software, 
and flow fields of plume impinging on different diversion troughs were compared and analyzed. It is found that, after 
LOX/kerosene engine exhaust plume vertically impacting on the diversion trough without deflector, plume which 
flows along the side surface of the diversion through and plume reflected upward would bring high temperature 
effect on the environment directly above the diversion trough. The wedge deflector would avoid the high 
temperature effect of plume on the environment directly above the diversion trough, while the conical deflector only 
could avoid the high temperature effect of plume reflected upward, there was still room for improving avoiding 
plume flowing along the side surface of the diversion trough. Compared with the plume impact on the diversion 
trough without deflector, the maximum pressure of the wedge deflector and the conical deflector respectively 
increase 10.31% and 33.81%, the plume impact maximum temperature of the wedge deflector and the conical 
deflector respectively decrease 4.04% and 8.95%. 
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0 引言 
液 氧 煤油 发 动机 因为 具有 推力 大 、 无 毒 、 无 污 


流 场 特性 仿真 研究 。Daniel C. Allgood 等 人 00 针 对 战 
神 系列 火箭 发 射 阶段 ， 对 氧 氧 发 动机 尾 焰 冲 击 导 流 
模 效 应 进行 了 CFD 仿真 计算 。Bruce T. Vu 等 人 0 


染 等 特点 已 经 成 为 液体 火箭 发 动机 的 发 展 方向 。 由 
于 火箭 在 发 射 过 程 中 喷射 出 高 速 、 高 温 的 尾 焰 燃气 ， 
在 发 射 初始 阶段 和 试验 平台 试车 过 程 中 发 动机 距离 
地 面 平台 较 近 ， 尾 焰 燃 气 的 热 冲击 和 动力 冲击 效应 
不 仅 可 能 引起 发 射 装 置 的 振动 响应 和 飞行 器 飞行 的 
初始 扰动 ， 甚 至 会 威胁 地 面 平台 上 其 他 设备 ("9， 因 


采用 CFD 计算 研究 了 航天 飞机 主 发 动机 和 固体 燃料 
助 推 器 尾 焰 对 有 忽 形 导 流 装置 导 流 槽 的 冲击 流 场 。 

本 文 提 出 了 一 种 液 氧 煤油 发 动机 尾 焰 冲击 导 流 
流 场 计算 方法 ， 首 先 计 算 发 动机 内 流 场 ， 然 后 以 


焰 冲 击 流 场 ， 该 方法 充分 考虑 了 燃烧 室内 燃烧 和 


模 
计算 获得 喷 管 喉 部 截面 参数 作为 入 口 边 界 条 件 计 算 
忆 
复 


此 有 必要 就 尾 焰 燃气 对 发 射 平 台 的 冲击 效应 进行 丰 


完 B。 在 航天 器 设计 过 程 中 ,通常 通过 计算 研究 火箭 
发 动机 尾 焰 对 喷 管 出 口 下 游 临 近 结 构 冲 击 的 热 影响 
3。 火箭 发 射 场 和 发 动机 实验 台 采 用 在 平台 下 面 建 
设 导 流 槽 的 方法 引导 尾 焰 流 向 ， 以 减缓 尾 焰 对 平台 
和 装置 设施 的 冲击 效应 ， 为 降低 导 流 权 的 冲击 烧 刨 
效应 、 提 高 导 流 槽 的 安全 性 和 使 用 寿命 ， 科 研 工作 
者 不 断 研究 不 同形 状 的 导 流 装置 。 

近年 来 ， 国 内 外 有 关 火 箭 发 动机 尾 焰 冲 击 导 流 
装置 的 数值 研究 不 断 深入 , 马艳丽 等 人 BI 针对 固体 火 
箭 发 动机 尾 焰 对 发 射 平台 冲击 进行 了 数值 研究 。 


燃 反 应 对 尾 焰 冲 击 的 影响 。 基 于 该 方法 采用 CFD 
软件 针对 液 氧 煤油 发 动机 尾 焰 冲击 效应 进行 三 维 数 
值 计 算 ， 研究 了 尾 焰 冲 击 不 同 结构 形状 导 流 装置 的 
流 场 特性 和 流动 规律 。 


1 物理 模型 和 计算 方法 


1.1 物理 模型 

本 文 以 液 氧 煤油 发 动机 为 研究 对 象 ， 推 进 剂 总 
流量 为 390 kg/s， 液 氧 /煤油 混合 比 为 2.4， 为 富 燃 燃 
烧 工 作 状 态 。 研 究 的 导 流 槽 及 导 流 装置 结构 如 图 1 
所 示 。 图 1 (a) 为 无 导 流 装置 的 简单 导 流 槽 ， 由 底 


Mayer 等 人 名 开发 程序 ， 计 算 研 究 了 基于 液体 火箭 发 
动机 尾 焰 冲击 的 传 热 特性 。Seiji Tsutsumi 等 人 中 针对 
固体 火箭 发 射 初始 阶段 冲击 流 场 开展 仿真 计算 ， 厂 
究 了 尾 焰 冲击 不 同形 状 导 流 装置 的 噪音 特性 。 
Michael 等 人 B9 开 发 固体 火箭 发 动机 冲击 流 场 仿真 
计算 程序 ， 通 过 建立 水 平 试验 台 ， 采 用 实验 方法 对 
开发 的 计算 程序 进行 验证 。Jeffrey 等 人 00 采 用 CFD 
软件 开展 火箭 尾 焰 冲击 有 圆锥 体 导 流 装置 导 流 槽 的 


面 、 侧 面 和 坡 面 组 成 。 发 动机 喷 管 出 口 直 径 为 DD， 
导 流 槽 结构 尺寸 为 侧面 高 度 a=7.64D (发 动机 喷 管 出 
口 与 导 流 槽 底面 距离 )， 底 面 横向 宽度 b=6.37D， 坡 
面 长 度 c=12.74D， 底 面 纵 向 宽度 4=6.37D。 图 1 (b) 
为 图 1 (a) 基础 上 在 导 流 槽 底面 建设 锦 形 导 流 装置 ， 
图 1 (c) 为 图 1 (a) 基础 上 在 导 流 槽 底面 建设 圆锥 
体 导 流 装 置 ， 锦 形 和 圆锥 体 曲 面 设 计 参 考 文献 [13]。 


人 而 合 面 。 发 动机 喷 管 0 
坡 而 Wp 坡 面 4 坡 面 
个 jm 7 se 一 圆锥 体 导 流 装置 
C 
业 b pe 


(a) 无 导 流 装 置 


(b) 有 铀 形 导 流 装 置 


(c) 有 圆锥 体 导 流 装置 


图 1 导 流 装置 和 导 流 槽 结构 


Fig.1 Structure size of diversion troughs and deflectors 


1.2 化 学 反应 机 理 


将 实际 煤油 燃烧 过 程 简 化 为 煤油 的 蔡 代 燃 料 C1i2H23 


文献 [1 入 对比 分析 了 液态 煤油 / 气 毛 / 气 氧 和 气态 
煤油 / 气 氧 / 气 氧 发 动机 内 部 流 场 计算 结果 ， 结 果 表 明 
两 种 计算 方法 都 可 以 获得 正确 流 场 计算 结果 。 本 文 
考虑 复 燃 对 尾 焰 冲 击 效应 的 影响 ， 采 用 简化 的 单 步 
总 包 反 应 


C,H;,+17.750, — >12CO,+11.5H,O (1) 


单 步 氧化 生成 完全 反应 产物 H2O 和 CO， 数 值 计算 
中 涉及 到 的 物质 有 煤油 蒸气 、 氧 气 、 二 氧化 碳 和 水 
蒸气 四 种 气体 。 采 用 单 步 总 包 化 学 反应 计算 得 到 的 
温度 相 比 实际 温度 会 偏 高 ， 但 是 计算 效率 却 非常 高 ， 
且 本 文 的 研究 目的 在 于 对 比分 析 不 同 导 流 槽 的 降 
温 性 能 ， 采 用 单 步 总 包 化 学 反应 可 以 得 到 正确 的 结 


论 。 
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1.3 控制 方程 与 求解 模型 
采用 多 组 分 化 学 反应 的 守恒 型 三 维 N-S 方程 作 


OU OF-F) 8(G-G) AH-H.) 
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为 模型 的 流动 、 物 质 与 能 量 交 换 以 及 燃烧 控制 方程 ， 
其 通用 形式 为 


At Ox Oy 


式 中 U 为 守恒 变量 向 量 ，1 为 时 间 变 量 ， 下 ， 
G ， 万 为 对 流 项 向 量 ，F 厂 ，G,， 顾 ,为 粘性 项 向 
量 ，J 为 源 项 向 量 。 上 式 分 别 为 质量 方程 、X，y， 
Zz 方向 动量 方程 、 能 量 方程 和 各 组 分 方程 。 
苗 述 火箭 发 动机 燃烧 流动 过 程 的 N-S 方程 是 非 
线性 强 相 互 而 合 的 偏 微分 方程 组 ， 目 前 压力 隐 式 算 
子 分 裂 算 法 (Pressure Implicit with Splitting of 
Operators，PISO ) 已 成 功 应 用 于 计算 火箭 发 动机 燃 
伐 流 动 过 程 05。 应 用 有 限 体积 法 离散 尾 焰 流 场 控制 
方程 ， 以 realizable Kk-é 双方 程 模型 和 有 限 速 率 / 涡 
耗 散 模型 求解 流 场 N-S 控制 方程 进行 计算 ， 使 用 
Arrhenius 公式 计算 化 学 源 项 。 
1.4 计算 网 格 与 边界 条 件 

基于 发 展 较为 成 熟 的 网 格 处 理 软件 对 网 格 进行 
处 理 , 由 于 发 动机 喷 注 器 面板 上 喷嘴 分 布 具有 1/3 对 
称 特性 ， 为 提高 计算 效率 ， 采 用 结构 化 网 格 的 方法 
处 理 本 文中 发 动机 内 部 流 场 三 维 计 算 区 域 ， 生 成 网 
格 及 边界 条 件 如 图 2 所 示 ; 由 于 尾 焰 冲击 导 流 槽 及 
导 流 装置 流 场 具 有 1/4 对 称 特性 ， 为 提高 计算 效率 ， 
选取 1/4 对 称 流 场 作为 尾 焰 冲 击 流 场 的 计算 区 域 , 生 
成 网 格 及 边界 条 件 如 图 3 所 示 。 经 过 网 格 独 立 性 验 
证 ， 最 终 选择 97 万 网 格 方案 和 111 万 网 格 方案 分 别 
进行 发 动机 内 流 场 和 尾 焰 冲击 流 场 计算 。 


图 2 发 动机 内 部 流 场 计 算 三 维 网 格 及 边界 条 件 
Fig.2 Three dimensional mesh and boundaries of engine internal 
flow field 
发 动机 内 部 流 场 计算 边界 条 件 具 体 定义 为 : 
mass flow inlet 边界 为 1.178 kg/s， 氧 /煤油 混合 比 为 


Oz = 
2.4; wall 边界 为 壁面 无 滑 移 ，pressure outlet 边界 为 
温度 300 K， 压 力 101325 Pa。 尾 焰 冲 击 流 场 计算 边 
界 条 件 具体 定义 为 : mass flow inlet 边界 通过 发 动机 
内 部 流 场 计算 确定 ; pressure far field 边界 为 温度 300 
K， 压 力 101325 Pa; wall 边界 为 壁面 无 滑 移 。 


pressure far field 


symmetry 
图 3 尾 焰 冲 击 流 场 计 算 三 维 网 格 及 边界 条 件 


Fig.3 Three dimensional mesh and boundaries of gasflow-guided 


channel and deflector 
2 计算 结果 与 分 析 


2.1 模型 验证 

使 用 本 文 模型 计算 文献 [16] 给 定 工 况 下 尾 焰 冲 
击 流 场 , 并 将 结果 与 之 进行 对 比 ,不 同 压力 比 (NPR= 
喷 管 总 压 / 喷 管 出 口 处 环境 压力 〉 工 况 下 计算 结果 与 
实验 测量 结果 对 比如 图 4 所 示 ， 左 边 是 本 文 计 算得 
到 流 场 速度 云图 ， 右 边 是 文献 中 纹 影 测量 流 场 结果 。 
喷 管 出 口 直径 4=2.54 cm， 喉 部 上 游 喷 管 收 缩 段 的 收 
缩 比 约 为 5， 地面 平板 与 喷 管 出 口 间 距离 h=24。 图 4 
(a) 是 NPR=2.5 时 流 场 对 比 图 , 流 场 中 存在 斜 激 波 ; 
图 4 (b) 和 图 4 (c) 分 别 是 NPR=3.7 和 NPR=5.0 
时 流 场 对 比 图 ， 流 场 不 再 产生 和 斜 激 波 而 产生 了 马赫 
盘 现象 ， 并 且 马 赫 盘 的 尺寸 随 着 相互 作用 强度 的 增 
大 而 增 大 ; NPR=3.7 时 , 不同 wd 工 况 下 冲击 平板 表 
面 压力 分 布 实验 对 仿真 对 比如 图 5 所 示 。 计 算得 到 
流 场 结构 、 平 板 表面 压力 分 布 与 实验 测 得 流 场 结构 
吻合 较 好 ， 证 明了 尾 焰 冲击 流 场 计算 模型 的 准确 性 。 


二 


201711.00618v1 


国 
国 


naXiv 


[1 


ChinaXiv 合 作 期 刊 


(a) NPR=2.5 (b) NPR=3.7 (ec) NPR=5.0 
图 4 不 同 压力 比 工 况 下 流 场 对 比 
Fig.4 Flow field contours of different NPR 
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(a) 实验 测量 结果 (b) 本 文 计 算 结 果 


图 5 冲击 平板 表面 压力 分 布 对 比 (NPR=3.7) 
Fig.5 Pressure distributions comparison between experiment and simulation on the plane surface (NPR=3.7) 

2.2 发 动机 内 部 流 场 计算 结果 与 分 析 以 发 动机 内 部 流 场 计 算得 到 喷 管 喉 部 截面 参数 
首先 ， 本 文 对 液 氧 煤油 发 动机 内 部 燃烧 流 场 进 作为 入 口 边界 ， 采 用 单 步 总 包 化 学 反应 描述 复 燃 反 
行 计算 ， 获 得 温度 和 速度 流 场 参数 分 布 如 图 6 所 示 ， ”应 过 程 ， 计 算 获得 了 液 氧 煤油 发 动机 尾 焰 冲 击 无 导 
氧 和 煤油 混合 富 燃 蒸气 在 发 动机 内 进行 充分 的 化 学 ，” 流 装置 、 有 狗 形 导 流 装置 、 有 圆锥 体 导 流 装 置 三 种 
反应 生成 二 氧化 碳 和 水 ， 喷 管 喉 部 和 扩张 段 流 场 参 “ 导 流 槽 的 三 维 流 场 ， 图 7 (a) -图 7 〈c) 分 别 为 尾 焰 
数 具 有 轴 对 称 特征 。 冲击 无 导 流 装置 导 流 槽 的 压力 、 温 度 和 速度 流 场 ， 
尾 焰 燃气 神 击 导 流 槽 底面 之 后 主要 治 着 导 流 模 的 坡 
面 、 导 流 槽 侧面 和 由 垂直 冲击 导致 的 向 上 反射 三 个 
方向 流动 ， 沿 着 导 流 槽 侧面 向 上 流动 和 向 上 反射 的 
燃气 都 将 对 导 流 槽 底面 正 上 方 环境 产 生 高 温 影响 。 
图 8 (a) -图 8 (c) 分 别 为 尾 焰 冲击 有 铅 形 导 流 装置 
ME 导 流 槽 的 压力 、 温 度 和 速度 流 场 ， 尾 焰 燃 气 直接 冲 

198385 击 钢 形 导 流 装置 顶端 之 后 主要 治 着 导 流 装置 引导 的 
坡 面 流动 ， 避 免 了 尾 焰 燃气 对 导 流 槽 底面 正 上 方 环 
境 的 高 温 影响 。 图 9 (a) -图 9《〈c) 分 别 为 尾 焰 冲 击 
有 圆锥 体 导 流 装 置 导 流 槽 的 压力 、 温 度 和 速度 流 场 ， 
尾 焰 燃气 直接 冲击 圆锥 体 导 流 半 置顶 端 之 后 沿 着 导 
流 装 置 引 导 方向 流动 ， 避 人 免 了 垂直 冲击 导 流 模 底面 
而 引起 的 燃气 向 上 反射 ， 但 由 于 圆锥 体 导 流 装 置 引 
导 一 部 分 燃气 治 着 导 流 槽 侧面 流动 ， 因 此 只 是 在 一 
定 程度 上 减缓 了 尾 焰 燃 气 对 导 流 槽 底面 正 上 方 环境 

(b) Velocity 的 高 温 影响 ， 还 具有 改进 的 余地 。 

6 发 动机 内 部 温度 与 速度 分 布 
Fig.6 Temperature and velocity contours in engine 


2.3 尾 焰 冲击 流 场 计算 结果 与 分 析 


六 


现 


201711.00618v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


m/s 
2744.78 
2058.58 
1372.39 
686.194 
0 


(a) 压力 (b) 温度 (c) 速度 
图 7 无 导 流 装 置 尾 焰 冲 击 流 场 


Fig.7 Impact flow field contours without deflector 
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(a) 压力 (b) 温度 (c) 速度 
图 8 有 锦 形 导 流 装置 尾 焰 冲击 流 场 


Fig.8 Impact flow field contours with the wedge deflector 
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(a) 压力 (b) 温度 (c) 速度 
图 9 有 圆锥 体 导 流 装 置 尾 焰 冲击 流 场 
Fig.9 Impact flow field contours with the conical deflector 

10 (a) -图 10〈c) 分 别 为 导 流 槽 底面 、 钙 形 ” 导致 发 动机 喷 管 出 口 与 冲击 表面 之 间距 离 缩短 了 ; 
导 流 装置 表面 和 圆锥 体 导 流 装 置 表面 上 压力 分 布 ， 二 是 导 流 装置 的 形状 导致 最 大 受 力 面 面积 减 小 了 ， 
锦 形 导 流 装置 和 圆锥 体 导 流 装置 表面 最 大 压力 相对 “并且 圆 锥 体 导 流 装 置 的 最 大 受 力 面 面 积 相对 忽 形 导 
无 导 流 装置 的 导 流 槽 底面 分 别 增 大 了 10.31% 和 ， 流 装 置 更 小 。 
3.81%， 这 是 因为 : 一 是 在 导 流 槽 底面 上 的 导 流 装置 


Y Y 
pressure: 84328.3 233242 382155 531069 679982 | pressure: 29137.5 209378 389618 569859 750099 Pa pressure: 34754.1 253539 472324 691109 909894Pa 
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(a) 无 导 流 装置 (b) 有 忽 形 导 流 装置 (c) 有 圆锥 体 导 流 装 置 
图 10 导 流 槽 底面 和 导 流 装置 表面 压力 分 布 
Fig.10 Pressure contours of the bottom surface and the deflector 
喷 管 喉 部 到 冲击 平面 间 中 心 轴线 上 参数 分 布 如 ”于 导 流 装置 能 够 很 好 的 引导 尾 焰 燃 气流 动 ， 相 对 尾 
图 11 所 示 ， 图 11 (a) -图 11(c) 分 别 为 压力 、 速 ”” 焰 燃气 冲击 无 导 流 装置 的 导 流 槽 ， 锦 形 导 流 装 置 和 
度 和 温度 分 布 ， 尾 焰 燃 气 冲击 到 导 流 槽 和 导 流 装置 。 圆锥 体 导 流 闭 置 表面 受 尾 焰 燃 气 冲击 的 最 高 温度 分 
表面 时 速度 又 降 为 0， 导 致 压力 和 温度 又 然 升 高 ; 由 ” 别 降低 了 4.04% 和 8.95%。 
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Fig.11 Parameters distribution on the axis 


3 结论 


通过 研究 得 到 如 下 结论 : 

(1) 液 氧 煤油 发 动机 尾 焰 燃气 垂直 冲击 无 导 流 
装置 的 导 流 槽 之 后 同时 沿 着 导 流 槽 的 坡 面 和 侧面 ; 
动 ， 沿 着 导 流 槽 侧面 流动 和 向 上 反射 的 燃气 会 对 
流 横 底 面 正 上 方 环 境 产生 高 温 影 响 。 

(2) 液 氧 煤油 发 动机 尾 焰 燃 气 冲击 有 锦 形 导 流 
装置 的 导 流 槽 之 后 沿 着 导 流 装置 和 导 流 覃 坡 面 流 
动 ， 避 免 了 沿 着 导 流 权 侧 面 流动 和 向 上 反射 的 燃气 
对 导 流 覃 底面 正 上 方 环境 的 高 温 影响 。 

(3) 液 氧 煤油 发 动机 尾 焰 燃 气 冲 击 有 圆锥 体 导 
流 装 置 的 导 流 槽 之 后 沿 着 导 流 装置 、 导 流 槽 坡 面 和 
侧面 流动 ， 一 定 程度 上 减缓 了 燃气 对 导 流 横 底 面 正 
上 方 环境 的 高 温 影响 ， 在 避免 燃气 沿 着 导 流 槽 侧 
流动 方面 仍 有 改进 余地 。 

(4) 相对 尾 焰 燃 气 冲 击 无 导 流 装置 的 导 流 模 ， 
由 于 导 流 装置 导致 发 动机 喷 管 出 口 与 导 流 冲击 表面 
之 间距 离 缩短 了 、 最 大 受 力 面 面 积 减 小 了 ， 锦 形 导 
流 装置 、 圆 锥 体 导 流 装 置 表面 最 大 压力 分 别 增 大 了 
10.31% 和 33.81%; 由 于 导 流 装置 能 够 很 好 的 引导 尾 
焰 燃 气流 动 ， 钢 形 导 流 装 置 、 圆 锥 体 导 流 装 置 受 尾 
焰 人 燃气 冲击 的 最 高 温度 分 别 降低 了 4.04% 和 8.95%。 
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